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Abstract 

A new method to deal with reduced dynamics of open systems by 
means of the Schrodinger equation is presented. It allows one to con- 
sider the reduced time evolution for correlated and uncorrelated initial 
conditions. 
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1 Introduction 



The study of time evolution of quantum open systems plays an important role 
in quantum information. The interaction between the system and the reservoir 
leads to phenomena of dissipation and decoherence [TJ [5] • These phenomena can 
be often modelled through the standard technique of quantum Markov processes 
in which the density matrix of an open system satisfies the Markovian master 
equation [SKUJS]. However, in the description of complex quantum systems one 
encounters a complicated non-Markovian behaviour [1,6-12]. Usually, non- 
Markovian evolution is studied under the assumption that the initial system- 
reservoir state is uncorrelated. 

The reduced dynamics can be considered in the cases when the state of 
the total system is mixed or pure. In the latter case, under the assumption of 
the existence of some symmetries one can solve the Schrodinger equation and 
deduce the reduced dynamics [13-21]. In the present paper a method to deal 
with reduced dynamics by means of the Schrodinger equation is presented. It 
allows one to study the reduced time evolution for correlated and uncorrelated 
initial conditions, and also makes it possible to consider models beyond the so 
called rotating wave approximation. 



2 A Class of Models 

We begin with the description of a class of models of open quantum systems for 
which the reduced dynamics can be studied in the case when memory effects as 
well as some initial correlations between the system and the reservoir are taken 
into account. The method presented here is an extended version of [22] . 

Let us consider a d-level system coupled to the boson field. The associated 
Hilbert space is C d <S> 0- where T is the bosonic Fock space. Let |1), . . . , \d) be 
a fixed orthonormal base in C d . The Hamiltonian of the system is assumed to 
have the form 

d 

H = lr + ld®H f + V (2.1) 

i=l 

with 

Hf = J dkui(k)a*(k)a(k) (2.2) 
while V is the interaction which is typically chosen in the form 

d 

V = K> 01 »<**(/«) +kc (2-3) 

Here a(k), a*(k), k € M. d are annihilation and creation operators on T and the 
following notation is used 

a(f) = [dk a(k)f(k) , (2.4) 
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where 

f(k)eL 2 (R d ). (2.5) 

It is clear that the Schrodinger equation with the above Hamiltonian cannot be 
solved since one has to deal with the whole Fock space. 

In order to restrict the reservoir degrees of freedom the interaction Hamil- 
tonian has to be simplified. Typically, one goes over to the so called V^-modcls 
[14], i.e. one assumes that not all transitions are allowed. More precisely one 
assumes that 

fij^O, i = 2,...,d, (2.6) 

and all the remaining form-factors vanish. In this case there are invariant sectors 
in <C d ®F. In each sector the corresponding Schrodinger equation can be solved. 
However, there is a different way to simplify the interaction Hamiltonian. Let 
|f2) be the vacuum state of the boson field, and let 

P = l s ®|fi)(fi| (2.7) 

be the projection operator on C d <S> T. Following [22] let us introduce the 
approximation V of V as follows 

V — ► V = PqV + VP - P VP . (2.8) 

In the case of the interaction V given in the form (|2.3p , one finds 

d 

V= \i)ti\®a*(fij)\tyM+h.c. (2.9) 
In this case the only possible transitions in T due to V are 

|n> «— o*(^-)|n). 

It should be pointed out that the above approximation makes 
beyond the rotating wave approximation. 



(2.10) 
it possible to go 



3 Time Evolution and Reduced Dynamics 

Let us consider a d-level system coupled to a bosonic reservoir for which the 
interaction Hamiltonian is given by (|2.9|) . If the state \<j>t) of the coupled system 
has the form 

d 

\4>t) - 52mci{t)\n) + a'(gi)\Sl)], (3.1) 
i=i 

the Schrodinger equation is reduced to a closed system of equations for Ci(t), . . . , Cd(t) 
and gl(k),...,gf(k), 

d 

Ck(t) = -ie k c k {t) - i};(fik, g\) (3.2) 

d 

g\ = -i( ei +u;(k))g l t -iJ2fin(k)c n (t), (3.3) 

n=l 
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where 

(j,g) = I dkf(k)g(k) (3.4) 



is the scalar product in L 2 (M. ). Moreover, the normalization condition 

(&|&) = (MM = 1, (3.5) 

takes the form 

d d 

Y^[h(t)\ 2 + (9l9l)} = E^(°)| 2 + ^,^)] = 1. (3.6) 

i=l i=l 

Equation (|3 . 3|) can also be written in the integral form 



ft 



(k) = e- u ^+^g^(k)-i / d S e-^- s )^ + ^( fc »E/ mn (fc)c„( S ). (3.7) 



71=0 



Inserting (|3.7p into (|3.2p one finds a closed system of equations for c\(t), . . ., 
Cfe(t) = -ie k c k (t)- J ds^M k i(t-s)ci(s)+G k {t), (3.8) 



o 

where 



;=i 



and 



Mki(t) = X) / dk f™kf™i e ~ lt{eMk)) ( 3 - 9 ) 

m=l ^ 

G fc (t) = -»X J dkf mk (k)g^(k)e- u ^ +u ^\ (3.10) 



The M k i(t) describe the reservoir correlation functions, while G k (t) corresponds 
to correlations between the reservoir and the initial state. Under the assumption 
that the initial state is given in the form 



d 

bo) = $>>® [ci(o)|n) + o*(flj)|n)] (3.ii) 



the total state \(j>t) is given by (|3.1[) in terms of the solution of the equations 
|, (|3.7p . It follows from (|3.1[) that the reduced evolution is given by 

d 

Trfl|0 t ><&| = X Ii>01fe(*)ci(*) + G^,ffj)]- (3-12) 
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In order to find explicitly the reduced evolution it is clear that the following 
correlation functions 



(t) = / dkf mn (k)f pq (k)e- iu W (3.13) 
W) = J dkf mn {k)g p {k)e-^ (3.14) 
c p , q {t) = [ dkg p {k)gl{k)e-^ (3.15) 



should be known. 

From the above definitions of the correlation functions, it follows that they 
satisfy the following positivity condition 

^ ^ ^ ^ ®"mn,pq{t'ot ^ ft)^mn (^a) %pq (^/?) 
ctrnn (3pq 

^2 ^2 C ^«(*a ~~ ^)^r (*a) (^) + (3-16) 

ar (3s 

(t a ) x r (tp) + c.c. > 

amn fir 

for all sequences of complex numbers x mn (ti), Xmnfo), x m n{t^), . . . and x p (ti), 
x p (t 2 ),.... 

By the Bochner theorem the correlation functions (|3.13p (|3.15[) can be pre- 
sented in the form 

oo 

a mn , pq {t) = [duje~ ZUJt J mn7 p q (uj) 7 (3-17) 
b mn , P {t) = I duje~ ZbJt J mn ,p(u) , (3.18) 



^ pq (t) = I dwe- iut J pq (u). (3.19) 



The matrices J m n,p q {w), Jmn,p(u), J pq (to) define the spectral density matrix 
such that for all sequences {x mn }, {x p } and u) G K the relation 

^ ^ Jmn,pq{^)x 7nn Xp q -\- ^ J rs {u)^jX r X s -\- ^ Jm^i^pi^^^X^^Xr -\- C.C. ^ 
ran p q rs mnr 

(3.20) 

holds. Thus, choosing the spectral density matrix J(w) one can model the time 
evolution of the reduced state. 

One can distinguish the special type of initial conditions, 

d 

l<M = (El J ) c ^°)) ®l n >- ( 3 - 21 ) 
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i.e., the simple factorized ones. In this case (|3.8p have the form 

d * 

Cfe(i) - -ie k c k (t)~Yl / dsM kl (t - s)ci(s) , (3.22) 
l=i { 

while the relation (|3 . 7|) reads 

\ d 

= -i dsJ2fmn(k)e-^- s ^ +u ^c n (s). (3.23) 
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Equations (|3.22p can be solved using the Laplace transform method, and their 
solution can be written in the form: 

d 

c k (t) = J2 L kl(t)ci(0)- (3.24) 

i=i 

On the other hand, from (|3.23[) it follows that 

d 

(9t,9t) = Rk m j m (t)c m {0)c n (0) . (3.25) 

m,n— 1 

It is worth noting that in the case considered the correlation functions (|3.13p are 
needed. Using (|3.12p . (13.24)1 and (|3.25p , one finds that the reduced evolution is 
given in the form 

jm in 

where 

Sjm,in (t) = Lj m (t)Li n (t) + Rjm,in (t) . (3.27) 

The relation (|3.26p can also be rewritten in the form 

Tr R \<j> t ) (<f> t \ = ^2 A t(\i) 01)^(0)^(0), (3.28) 

ijm 

where the relation 

At (|i> 01) = H H Sn ^m (t) (3.29) 

inn 

defines a map 

A t : M d — » M d , t > , (3.30) 

which is completely positive and trace preserving by construction. 

The discussion of time evolution under initial conditions different from (|3 . 2 1 [) 
is, in general, very complicated and will be presented on a simple model. 
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4 Spin-Boson Model beyond the Rotating Wave 
Approximation 

The model is specified by the Hamiltonian 

H = lo<t+(t-+ ( dkuj{k)a*(k)a{k) + a-®a*(f)\n)(rL\+<r+®a*(h)\n)(ri\+h.c. 



(4-1) 

on C 2 ® J and can be obtained from the spin-boson Hamiltonian containing 
anti-resonance interaction using the method presented in Sect. 3. The wave 
function can be chosen in the form 

\cj> t ) = (c (i)|0) + Cl (i)|l)) ® \Sl) + |0> ® a*( 5t °)|fi) + |1}) ® a*^ 1 )^) (4.2) 



since the interaction in (14. jj allows only the transition between the vacuum state 
and 1-boson states. The Schrodinger equation is equivalent to the following set 
of equations for Co(i), c\(t) and g°(k), g\{k): 

co = -(h,gl), (4.3) 

ci = (*)-*(/, (4.4) 

S°(fc) = -iw(fc) St (fc)-t/(fc)ci(t)), (4.5) 

= -i(u + u(k))gi{k)+h{k)co(t)). (4.6) 

Equations (|4.4|) and (|4.5|) are the same as in the case of spin-boson model in 
the rotating wave approximation. On the other hand (|4.3D . (|4.6|l describe the 
effect of anti-resonance interaction. It follows from l|4.3[) - (|4.6[) that our system 
has two constants of motion: 

\co(t)\ 2 + (gt,gl) = \co(0)\ 2 + (sb>9o) = Po > (4.7) 

and 

\ci(t)\ 2 + (g° t ,g°) = |ci(0)| 2 + ( ff0 °, ffo °) = pa > (4.8) 
Moreover, the normalization condition 

<&|&> = (^ol^o) = 1, (4.9) 

takes the form 

Po+Pi = 1. (4.10) 
The reduced state Tr^^l^t) can be written as follows 

l 

p(t) = TtRfalfr) = ^2(i\j) Pij (t), (4.11) 



i,j=0 



where 



poo® = |c (£)| 2 + (.g0, fft °) (4.12) 
pn(i) = \ Cl (t)\ 2 + (gl,gt) (4.13) 

An(*) = M*) = co(*)^(*) + (fft,fft)- ( 4 - 14 ) 
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Using (|4.7[) and (|4.8[) the diagonal elements of p(t) can be written in the form 

poo(*) = |coW| 2 + |c 1 (0)| 2 -|c 1 (t)| 2 + ( ffo , 5o ) (4.15) 
p u (t) = l-poo(*). (4-16) 

It should be pointed out that the equations (|4.3p . (|4"1))) and (|4.4p . (|43| are 
independent. However, there exists coupling in terms of correlation functions 
which are needed to calculate the reduced evolution. Eliminating g®{k) and 
g\ (k) from (|4.5|) and (|4.6p the following equations for co(t) and c±(t) can be 
derived 

t 

c (i) = J ds m (t - s)co(s) + n (t) , (4-17) 
o 

t 

ci(i) = -iwci(t) - \ dsmx{t- s)c x {s) + m(t) , (4-18) 



where 



(i 



mo(t) = J dk \h(k)\ 2 e- it( - UJ+ ^ k) \ (4.19) 
m x (t) = / eft |/(/c)| 2 e- lMfc) (4.20) 



and 



no(t) = -i J dk h{k)gl{k)e- it{ui+w{ - k)) , (4.21) 
ni(t) = -i [ dk f(k)g° (k)e- u ^ k K (4.22) 



In order to find out an explicit form of p(t) a detailed knowledge of all the 
correlation functions is required. However, general properties of poo(i) and 
pn(t) can be deduced in the following manner. Suppose that Ao(t)co(0) and 
Ai(t)ci(0) are general solutions of homogeneous equations (|4.17|) and (|4.18|) . 
respectively. Then one has 

t 

Ca(t) = A a (t)c a (0) + J ds A a {t - s)n a (s) (4.23) 
o 

with a = 0,1 and 

Poo(t) - IcoWl 2 + |ci(0)| 2 - | Cl (i)| 2 + (glg°) . (4.24) 

It follows from (|4.23|) that the time dependence of pao(t) is related to the initial 
conditions. Moreover, the asymptotic behaviour of poo(t) can be easily found as 

lim p Q0 (t) - |ci(0)| 2 + (.9 ,5o ) (4-25) 

t — *oo 
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since, typically, one has 



lim c a (t) — , a = 0, 1 . 



(4.26) 



The formula (|4.25|) shows that the asymptotic state is determined by the ini- 
tial conditions, i.e., the equilibrium state does not exist. In the case g^k) = 
9o(k) = the time evolution is given in terms of a completely positive and trace 
preserving but not relaxing map. 

The spin-boson model is recovered if one puts h(k) = = g^k) = g\(t). In 
this case we have 
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